Indications of a Metallic Antiferromagnetic Phase in the 2 Dimensional U — t — t' Model 
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We present mean-field and quantum Monte Carlo results that suggest the existence of an itinerant 
antiferromagnetic ground state in the half-filled U — t — t' model in two dimensions. In particular, 
working at t'/t = —0.2 we found that antiferromagnetic long range order develops at U ci /t ~ 
2.5 ± 0.5, while a study of the density of states N(ui) and the response to an external magnetic field 
indicates that the system becomes insulating at a larger coupling 4 < U C2 /t < 6. 

PACS numbers: 71.30.+h, 75.10.Lp, 75.50.Ee 



The interest in metal- insulator transitions started sev- 
eral decades ago with the observation that nickel oxide 
(NiO), a transparent non-metal, should be metallic ac- 
cording to its electronic band structure. [jjj Afterwards, 
several models were proposed to study metal-insulator 
transitions, including the well-known Hubbard model. 
Brinkman and Rice found that when the on-site Coulomb 
repulsion U reaches U c = 1A5W, where W is the elec- 
tronic bandwidth, || the ground state should change 
from metallic to insulating. However, a variety of nu- 
merical and analytical studies have convincingly shown 
that the half-filled Hubbard model with electronic hop- 
ping between nearest neighbor sites has an insulating an- 
tiferromagnetic ground state for any finite value of the 
coupling U. The addition of a hopping along the plaque- 
tte diagonals to the Hubbard Hamiltonian destroys the 
nesting, and the possibility of a finite critical coupling 
(U c ) at half-filling is recovered. U c depends on the diago- 
nal hopping t' , and in previous work it has been obtained 
by monitoring the value of U at which antiferromagnetic 
long range order (AFLRO) develops. || Mean field and 
quantum Monte Carlo methods have been the techniques 
most used for this purpose. It is generally accepted that 
the half-filled U — t — t' model has a metallic ground state 
for U < U c and an antiferromagnetic insulating (AFI) 
state for U > U c . 

The aim of our work is the study of the metal-insulator 
transition (MIT) in the U — t — t' model using updated 
numerical techniques. We will monitor the development 
of AFLRO from the behavior of the spin correlations, 
and we will search for the onset of an insulating phase 
by studying the density of states and the response of the 
system to magnetic fields. The main result is that we 
have found indications of the existence of an intermedi- 
ate phase between the paramagnetic metal and the AF 
insulator. The phase can be characterized as an AF metal 
(AFM). The existence of two dimensional systems with 
non-ordered metallic, AFI and AFM phases has recently 
been discussed. Hi 



The U — t — t' Hamiltonian is given by 
H=-t ^ (c^c^ + h.c.)-!' Y, ( c Lc„, ff +h.c.) 
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where Cj creates an electron at site i with spin projec- 
tion a, nio- is the number operator, the sum (ij) runs over 
pairs of nearest neighbor lattice sites, and the sum (in) 
runs over pairs of lattice sites along the plaquette diago- 
nals. U is the on-site Coulombic repulsion, t the nearest 
neighbor hopping amplitude, t' the hopping amplitude 
along the plaquette diagonals, and /i is the chemical po- 
tential. In the following t — 1 will be used. 

As a first step, the Hamiltonian will be studied using 
the spin density wave (SDW) mean-field (MF) approxi- 
mation. |y Proposing as an Ansatz an AF ground state 
and following the standard Bogoliuvov procedure to di- 
agonalize exactly the resulting MF Hamiltonian we found 
two energy bands given by 
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e£ = —At' ' cosk x cosk y , 
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A is the MF parameter that, when finite, indicates that 
the ground state has AFLRO. Both A and n are obtained 
by solving the self-consistent equations: 
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FIG. 1 . SDW mean-field phase diagram of the U-t-t'- if' 
model at half-filling and T = 0, for t = 1 and t' = -0.2. 
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where /(#) is the Fermi function given by e pl +1 ■ 
the number of sites in the lattice and 
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n(k) = ±(1 - ^)/(£ k ) + ^(1 + =)/(2£). (8) 
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Eqs.(6-8) can be solved at any temperature T but since 
we are trying to study ground state properties we will 
work at T = 0. To study the MIT discussed before 
(n) = 1 is fixed in our calculations. When if = it is 
found that as soon as A becomes different from zero, in- 
dicating the existence of an AF ground state, a gap of size 
2 A opens between the two bands given by Eq.(2), and 
the chemical potential lies inside the gap. However, when 
t' ^ the shape of the bands is distorted, and even when 
AF has developed, the two bands overlap if A < 2\t'\. 
This occurs because the actual gap is defined by the sep- 
aration between the highest state in the lower band and 
the lowest state in the upper band. For t' < these 
states have k = (it/2,it/2) and k = (ir,0), respectively 
(for t' > the momenta are reversed). For example, 
for t' < 0, the lowest state in the upper band has energy 
A — A\t'\ while the highest state in the lower band has en- 
ergy —A. An analogous condition is obtained for t' > 0. 
Then, an effective finite gap exists if A > 2\t'\. There- 
fore, in the region where the bands overlap the ground 
state is antiferromagnetic, but it has metallic properties 
because the chemical potential cuts both the lower and 
upper bands. Unfortunately, solving the MF equations 
we found that the AFM is not stable in this case since A 
changes discontinuously from zero to 2\t'\ at a particular 
value of U which depends on t' . This is due to the fact 




(0,0) (Hit) (0,Jt) (0,0) (0,0) (Jljt) (0,11) (0,0) 

FIG. 2. SDW mean-field band structure of the half-filled 
U -t-t' -t" model for t = 1, t' = -0.2 and t" = 0.2, a) in 
the AFM phase with U = 3.57 and A = 0.8 and b) in the AFI 
phase with U = 5.1 and A = 2.0. The dashed line indicates 
the position of the chemical potential. 

that the energy as a function of A, in this case, has a 
two minima structure. Then, the AFM phase exists only 
at one point in parameter space. However, it is reason- 
able to expect that the AFM phase might be stabilized 
by including the effect of fluctuations beyond the MF 
approximation. Rather than adding these fluctuations, 
we intuitively believe that their effect could be mimicked 
by introducing longer range hopping terms in the kinetic 
energy. Thus, we have added to the MF calculation an 
additional hopping term between second nearest neigh- 
bors in the x and y directions with strength regulated by 
the parameter t". The MF equations are modified simply 
by replacing in Eq.(3) by 



-At' ' cosk x cosk y — 2t" (cos2k x + cos2k y 
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In this case we found that the overlap of the bands occurs 
if < A < 2\t'\+At". The introduction oft" stabilizes the 
AFM phase in the MF approximation as can be seen in 
the phase diagram presented in Fig.l for t' = —0.2. The 
circles indicate the values of U Cl where the ground state 
changes from paramagnetic metal to AF metal, i.e. where 
A becomes non-zero but the bands still overlap. The 
squares represent a second critical coupling, U C2 , at which 
the ground state becomes an AF insulator. This occurs 
when A becomes larger than 2\t'\ + At" and the bands 
no longer overlap. In Fig. 2. a the energy bands along 
certain directions in momentum space in the AFM phase 
are shown. As it was discussed above, it is clear that 
the two bands overlap and the position of the chemical 
potential, denoted with a dashed line, indicates that the 
top of the lower band is empty while the bottom of the 
upper band is filled. The energy bands in the AFI phase 
are presented in Fig.2.b. In this case, /i lies in the middle 
of the gap and the system is clearly an insulator. Then, 
at the MF level, the addition of longer range hopping 
terms stabilizes the AFM ground state. || 

To find further support for the existence of the AFM 
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FIG. 3. QMC spin-spin correlation C(r) = (S-S? +r )(-l) |r| 
for several values of U/t, T = t/6 on an 8x8 lattice at 
half- filling for t' /t = —0.2. Points without error bars have 
errors smaller than the size of the dots. 

ground state, Quantum Monte Carlo (QMC) techniques 
have been here applied to the study of the U—t—t' model, 
working at t' = —0.2 on 8 x 8 lattices and temperatures 
T = 1/8 and 1/6. For these values of the parameters 
the MF approximation (with t" = 0) predicts U c = 2.1. 
Previous Quantum Monte Carlo results suggested that 
AFLRO develops for U C1 « 2.5. j| We have indepen- 
dently analyzed the spin-spin correlation function C(r) 
as a function of distance for several values of U/t on an 
8x8 cluster. The results are shown in Fig. 3. Note 
that a finite tail is already developed for U/t — 3 but 
it is not present for U/t = 2. Thus, a critical coupling 
U c « 2.5 ± 0.5 is here estimated. ffl] Now let us focus 
our attention to the coupling [7 = 4 where AFLRO is 
clearly developed (Fig. 3). The next issue is whether the 
system is metallic or insulating for [7 = 4. To investigate 
this important point, the density of states N(lo) was cal- 
culated using the maximum entropy technique. S The 
stability of the results was checked by making four inde- 
pendent long runs for each set of parameters. In Fig. 4. a 
the results for U = 4, t' = —0.2 and = 6 are presented. 
Three main peaks are observed in the results: those to 
the left and to the right of n can be identified with the 
lower and upper Hubbard bands, while the peak at /u, cor- 
responds to the quasiparticle weight indicating that the 
system is metallic. This metallic behavior disappears as 
U increases and an effective gap develops. In Fig.4.b 
N(lj) for U = 6, t' = -0.2 and = 4 is shown. Here 
the finite temperature precursor of an insulating gap in 
the density of states has developed and the chemical po- 
tential lies inside the gap. Notice that these results are 
in agreement with those we obtained in Rcf. JoJ where 
the gap was studied by monitoring the behavior of the 
density (n) versus /i, as well as the spectral functions 
A(k, u). In Fig.4.c the mean field result for N{ui) in the 
AFM region is presented. It is remarkable that the three 
main peak structure is qualitatively similar to that found 
numerically for U = 4. The MF density of states in the 
AFI phase is presented in Fig.4.d where there is also good 



FIG. 4. QMC density of states N(u) for the U-t-t' model 
at half-filling for t' = -0.2 and a) U = 4, (3 = 6, on an 8 x 8 
lattice: b) same as (a) for (7 = 6 and j3 = 4; c) mean-field 
results on a 200 x 200 lattice at T=0 using t" = 0.2, A = 0.6 
and [i = —0.382; d) same as (c) with A = 1.8 and fj, — 0. The 
MF transition between AFM and AFI occurs for A = 1.2 

qualitative agreement with the QMC numerical results 
for U = 6. 

In spite of the good agreement between the SDW-MF 
and QMC numerical data suggesting the existence of an 
AFM phase and a metal-insulator transition somewhere 
in the interval 4 < U C2 < 6, we decided to study the lower 
limit for U C2 using an alternative technique. It is well 
known that when a metallic system is close to its critical 
coupling towards an insulator, a magnetic field can in- 
duce a metal-insulator transition. |10|Jll] ] The signature 
of the transition at low temperature is a discontinuity 
( "metamagnetic transition") in the magnetization m as 
a function of the magnetic field h for U slightly smaller 
than U C2 . For a coupling close to but smaller than U C2 
the system is still metallic at zero magnetic field. The 
magnetization curve follows a metallic behavior until the 
field reaches a critical value h c that drives the system to 
become an insulator and, thus, producing the metamag- 
netic transition. Since we are working at small but finite 
temperature, we do not expect to find a perfect disconti- 
nuity but instead a very rapid crossover. In fact, in Ref. 
]IT[ , where the Hubbard model in infinite dimension was 
studied, it was shown that only for T < 0.01 a proper 
first order phase transition is observed. However, defin- 
ing x = dm / dh as the magnetic susceptibility, at low but 
finite temperatures and at low magnetic fields we expect 
dx/dh > to be the signature of the metamagnetic tran- 
sition that should occur at lower temperature. On the 
other hand, dx/dh < would indicate normal behavior. 
pl[ For a standard Fermi liquid, m vs h should be linear 
for small magnetic fields (i.e., dx/dh = 0). With increas- 
ing h, the slope will decrease as m saturates to 1. The 
numerical calculation of m as a function of h is very dif- 
ficult because the behavior that we want to study occurs 
at very small magnetic field, i.e. h < 0.1. We found that 
the anomalous behavior is not observed for temperatures 
higher than T = 0.125. To obtain numerical values 
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FIG. 5. Magnetization m versus magnetic field h for the 
half-filled U -t-t' model, t' = -0.2 and U = 2, 3 and 4 at 
T=t/8. In the inset m versus h is shown for U — 4 and t' = 
(filled diamonds) and -0.2 (open diamonds). 

with small enough error bars we needed to perform about 
100,000 measuring QMC sweeps per point. In Fig. 5, m 
vs h for t' = —0.2 is shown at different values of U for 
= 8 on an 8 x 8 lattice, dx/dh < is observed for 
U = 2 and 3. However, for U = 4 dx/dh > is found. 
This suggests that at t' = —0.2, U = 4 is a lower bound 
for U C2 . It could be argued that a quantum antiferro- 
magnet will have dx/dh > because its magnetization is 
given by ui = x.h+ sgn h h 2 /Aire 2 |Q. However, for the 
small values of h used here the linear term prevails as can 
be seen in the inset of Fig. 5 where the magnetization as 
a function of the magnetic field is shown for U = 4 and 
t' = 0. Since antiferromagnetism is reduced by the addi- 
tion of a finite t' the upturn observed for t' = —0.2 can 
not be caused by staggered spin correlations. Then, by 
comparing the two curves in the inset of Fig. 5, it is clear 
that dx/dh becomes larger than zero due to the effect 
of the diagonal hopping. It appears that in the half- 
filled U -t-t' model when f' = -0.2, AFLRO develops 
at U Cl — 2.5 ± 0.5 but the system becomes an insula- 
tor for U C2 larger than 4. These numerical results again 
suggest that the AFM region exists, and it is actually 
broader than what was predicted at the mean field level. 
Notice that the existence of AFM phases is not just an 
academic curiosity. Magnetically ordered metallic phases 
have been observed experimentally. E| In particular, an 
AFM phase was found in organic k - (BEDT - TTF) 2 X. 


Summarizing, in this paper we have provided numer- 
ical and analytical results that suggest the existence of 
an antiferromagnetic metallic ground state in the two di- 
mensional U — t — t' model. 
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